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ABSTRACT 



A sequence of pure-hydrogen C05BOLD 3D model atmospheres of DA white dwarfs is presented for a surface gravity of log g = 
8 and effective temperatures from 6000 to 13,000 K. We show that convective properties, such as flow velocities, characteristic 
granulation size and intensity contrast of the granulation patterns, change significantly over this range. We demonstrate that these 
3D simulations are not sensitive to numerical parameters unlike the ID structures that considerably depend on the mixing-length 
parameters. We conclude that 3D spectra can be used directly in the spectroscopic analyses of DA white dwarfs. We confirm the result 
of an earlier preliminary study that 3D model spectra provide a much better characterization of the mass distribution of white dwarfs 
and that shortcomings of the ID mixing-length theory are responsible for the spurious high-log g determinations of cool white dwarfs. 
In particular, the ID theory is unable to account for the cooling effect of the convective overshoot in the upper atmospheres. 

Key words, convection — hydrodynamics — line: profiles — stars: atmospheres — white dwarfs 



$h \ 1. Introduction 



The multidimensional radiation-hydrodynamics (RHD) treat- 
ment of convective motions in hydrogen-atmosphere DA white 
dwarfs was pio neered by the Kiel group some fifteen years ag o 

\ dLudwig et alj 11994 ISteffen et al.l 1 1 995b iFrevtag et alj [l996). 
These studies have revealed that the differences between spec- 

\ tra computed on the basis of 2D RHD simulations and 
those computed from standard ID hydrostatic models were 
small. However, recent research developments make it advis- 
able to pursue the computation of RHD models. First of all, 
three-dimensional RHD computations (hereafter 3D simula- 
tions) became accurate enough to be the reference in quan- 
titative spectroscopy, in particular for the abundance determi - 
nations in the Sun dAsplund et all 120091: ICaffau et al] 1201 ll) . 
Secondly, the spectroscopic technique in DA white dwarfs, 
which consists in comparing the observed line profiles of the 
B aimer series with the predictions of detailed model atmo - 
spheres (Weidema nn & Koested fl980: Berger onet al.1 1992ah 
becam e more widely used, with larger samples (lEisenstein et alj 
2006 ) and high er signal-to-noise ratio (S/N) observations 
(Gia nninas et al.l 1201 lh . For high quality observations (S/N > 
50), it is often quoted that the relative uncertainties in the atmo- 
spheric parameters, the effective temperature (T e tf ) and surface 
gravity (log g, in cgs units), are of the order of 1% dLiebert et al.l 
120051: lKoesteretal]|2009bl) . As a consequence, white dwarfs 
have been used as spectroscopic standards t o achieve hi gh- 
precision calibrations at the 1% level (see, e.g.. lBohiinll2000l) . 

It has been shown long ago that the spectroscopically de- 
termined masses of cool DA white dwarfs (r e ff < 13,000 K) 
were as much as 20% higher than those of hotter DA white 
dwarfs dBergeron et al.|[l990l Il992al) . This discrepancy is now 
observed in the spectroscopic mass distribution of all surveys 



of DA white dwarfs dTremblav et al.l 1201 lat iGianninas et al.l 
12011). Multiple attempts were made to ex plain this behaviour 
(Koe ster et alJl2009atlTremblav et al.l l2010) but were unsuccess- 
ful. However, it was suggested that because convection becomes 
significant in the photosphere almost exactly as the so-called 
high-log g problem become s apparent, the ID mixing-length 
theory dBohm-Vitensdl 19581 hereafter M LT) was likely the cul- 
prit. Therefore, iTremblav et al.l (1201 lbl hereafter Paper I) com- 
puted the first four non-grey 3D model atmospheres of DA white 
dwarfs (12,800 > T eS (K) > 11,300 and log g = 8) using the 
C0 5 BOLD RHD code dFrevtag et alJl2012h . The code is usually 
applied to stars with deep convective envelopes where it is not 
possible for the simulation domain to span the full convection 
zone. However, the convection zone in these hot white dwarfs is 
in fact thinner than the simulation domain in the vertical direc- 
tion. 

The aim of Paper I was to predict quantitatively how the at- 
mospheric parameters are changed by using 3D spectra instead 
of the standard ID spectra. As a consequence, the methodologi- 
cal approach was to do a differential analysis, by comparing both 
the spectra comp uted from 3D structures , and ID hydrostatic 
(LHD) structures (Caffa u& Ludwigll2007l) relying the same mi- 
crophysics and opacity bins as the 3D simulations , with standard 
ID spectra of white dwarfs (ITremblav et al.l201 lab . It was found 
that the 3D log g corrections had the correct amplitude to solve 
the high-log g problem for DA white dwarfs in the range of T e s 
studied in the paper. It was concluded that a weakness in the ID 
MLT theory was creating the high-log g problem. 

Several aspects of the 3D calculations can be improved in 
order to understand to which extend these simulations can solve 
the high-log g problem, and ultimately, to use them in the spec- 
troscopic analyses of white dwarfs. First of all, it is desirable 
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to compute cooler 3D simulations, since the high-log g problem 
is most prominent around T e s = 10,000 K (see, e.g., Fig. 1 of 
Paper I). Furthermore, it is necessary to have 3D simulations 
containing realistic microphysical properties in order to compare 
them directly to white dwarf observations. For instance, the four 
3D simulations computed in Paper I rely on an equation-of-state 
(EOS) and opacity table (including an opacity binning proce- 
dure) that have different microphysics in comparison to standard 
ID model atmospheres. 

We tackle these shortcomings in this work by producing an 
improved sequence of 12 DA model atmospheres, covering the 
range 13,000 > r eff (K) > 6000 at logg = 8. In Sect. 2, we 
discuss about the precision of our 3D simulations. In Sect. 3, 
we present the properties of our sequence of 3D computations. 
Model spectra are computed in Sect. 4 and improved 3D log g 
corrections are derived in relation to the high-log g problem. The 
conclusion follows in Sect. 5. 

2. Establishing absolute properties 

The main objective of our study is to produce 3D model at- 
mospheres that represent real stellar conditions as closely as 
possible with the aim of performing direct comparisons of 
these models with observations. We must beforehand under- 
stand the sources of the differences between these 3D struc- 
tures and standard ID structures, both computed with estab- 
lished well-known codes. In other words, we must carefully 
compare the codes producing these model atmospheres. Several 
groups have in the past computed grid s of ID model at- 
mospheres for DA white dwarf's (e.g.. iBergeron et al.l Il992al 
iFinlev et all Il997t iKowalski & Saumonl 120061) . Recent theo- 
retical developments include the improved modeli ng of the 
Lyman quasi-molecular satellites (lAllard et al.l 120041) . the cal- 
culation of the Ly-q red w ing opacity due to H2-H collisions 
(IKowalski & Saumonl |2006[) and the publication of new Stark 
broadening profiles (Tre mblav & Bergeronlt2009|) including in a 
consis tent way the non-ideal gas effects of iHummer & Mihalas 
( 1988). All of these improvements are now included in ou r ID 
model atmosphere code described in Tremblav et Iii1(l201 lal) and 
references therein. We note that the D. Koester model atmo- 
sphere code in cludes very similar physics at the time of publica- 



tion (s ee, e.g., Girven et al.l20lil) . Hence, while we specifically 



use the lTremblav et al.l(l201 lal) grid in this work, it is appropriate 
to describe them from now on as standard ID models. 

The treatment of convection is obviously a genuine differ- 
ence between the 3D CO s BOLD code and the standard ID 
codes. As a consequence, it is difficult to isolate other possible 
differences between both codes. This hurdle is removed by using 
the ID LHD code, in which microphysics and radiative trans- 
fer numerical schemes (e.g., the opacity binning) are identical 
to those included in the 3D simulations. In Fig. Q~l we compare 
our latest ID LHD structures and standard ID structures, which 
shows that we were successful in explaining and eliminating as 
much as possible the differences between both codes. Our im- 
provements made to the LHD code are described in this section. 

2.1. input microphysics 

In Paper I, the CO 5 BO LD and LHD struc tures were based on 
different microphysics dLudwig et al.l!994b in comparison to the 
one used in prevalent ID white dwarf codes. The former codes 
rely on pre-computed EOS and opacity tables as input, and it 
is fairly straightforward to replace these tables with ones that 
are computed with the same microphysics as in standard ID 



ID LHD (8 opacity bins) 

ID standard model (1813 frequencies) 
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Fig. 1 : Atmospheric temperature as a function of the logarithm of 
the Rosseland optical depth for DA models with effective tem- 
peratures of 8000 (bottom) and 12,000 K (top ). We show the 
standa rd ID structures (solid black line) from iTremblav et al.l 
d201 lal) . recomputed with a slightly higher numerical precision, 
and the improved ID LHD structures (dashed blue line) de- 
scribed in the text. 



models. This is the first step of our improvement of the LHD 
structures. We found that the non-ideal EOS typically used in 
white dwarfs has a negligible effect on the predicted structures 
(for r e ff > 6000 K and log tr < 3) compared to the ideal EOS 
utilised in Paper I. In the case of the improved opacities, the ef- 
fect is mild, and it is shown in Sect. 4.3 that the corresponding 
log g correction is of the order of 0.03 dex. 

We now look at the opacity binning approach, which is a 
known difference between the LHD and standard white dwarf 
codes. In Paper I, we used 7 band-averaged opacities to de- 
scribe the band-int egrated radiative t ransfer, based on the pro- 
cedure laid out in Nordlund (1982); Ludwig et aD (11994J) and 
Vogl er et al.l (12001 7 These include two bins dedicated to the 
Lyman pseudo-m olecular satellites (Alla rd & Kielkopi] [l982: 
All ar d et al]|2004l) . The ID models calculated in Tremblav et al.l 
( 12011 ah . "on the other hand, employ 1813 carefully chosen fre- 
quencies for the radiative transfer (which is sufficient for the 
broad lines in DA stars). 

Fig- H] presents the characteristic forming region (r,i = 1) 
of the emergent flux as a function of the wavelength. We rely 
on the Rosseland optical depth (tr) as the reference scale. 
Unsurprisingly, the lines are formed higher in the atmosphere 
than the continuum, since the opacity is larger at these wave- 
lengths. We see that the atmospheric region within -5 < 
log tr < is relevant for the line formation. 

We sort the wavelength-dependent opacities based on the 
Rosseland optical depth at which = 1. In other words, the 
boundaries of the opacity bins could be represented by hor- 
izontal lines in Fig. [2] We use a total of 8 opacity bins in 
our calculations, by applying thresholds in log tr given by 
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[oo, 0.0, -0.5, -1.0, -2.0, -3.0, -4.0, -oo] and adding one bin 
for the Lyman satellites. We find that it is necessary to use at 
least one bin for each dex in tr to reproduce the standard ID 
calculations in that tr region. We also split in two bins the criti- 
cal -1 < log tr < region to be on the safe side. In all but two 
bins a switching between Rosseland and Planck averages is per- 
formed at a band-averaged Rosseland optical depth of 0.35. In 
the two bins gathering the largest line opacities, the Rosseland 
mean opacity is used throughout. W e treat scattering as tru e ab- 
sorption. We have verified with the Tremb lav et al.l (1201 lab at- 
mosphere code that this approximation has a negligible effect on 
the structure of convective DA white dwarfs. Fig. [TJ shows that 
we can reproduce very well the established ID structures with 
our opacity tables. This is an important result, since it was not 
obvious that we could reproduce with accuracy the structures of 
standard ID models with only a few opacity bins. We note, in 
comparison, that heating rates in the Sun are wel l reproduced 
with a 9-bin/12-group scheme (Freytag et al. 2012|). 



ail 
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Fig. 2: Value of the Rosseland optical depth at which the plasma 
becomes optically thin (t,i = 1) at a particular wavelength for ID 
model atmospheres at 6000 (red, solid), 10,000 (black, dotted) 
and 12,000 K (blue, dashed). 



2.2. Treatment of convection 

One aspect that could explain differences between the LHD 
and standard structures is a difference in the physical parame- 
ters that are allowed to vary. Unfortunately, we discovered such 
incompatibility which is related to the MLT convection input 
parameters used to compute LHD structures. The LHD code 
setup in this study relies on the ML2/o- = 0.8 parameterisation 
of th e MLT that is used in the most recent white dwarf mod- 
els (Tremblav et al. 2010), where the calibration of the mixing- 
length parameter a results from a comparison of near-UV and 
optically determined T e ff. For completeness, we summarise in 
Table 1 the a,b and c parameters that are use d to define the 
ML2 treatment of the MLT equations (see , e.g., Mihalas 1978; 
Tass oul et all 1 19901: iBergeron etaT1ll992bl iKoester et alii 1994 
Ludwig et al.l [l999). There is a fourth parameter that was intro- 



duced early in the usage of the MLT theory (IVitenselll953l) to 
allow for optically thin convective cells (see eq. 7.72 of Mihalas 
1978). This parameter is defined as d in the white dwarf field 
(see eq. 4 of Berger on et al.l 19 92b). and a U paramete r, specified 



in a different way but recoveri ng the Mihalas! ( 19781) relation if 
fa = 2, is defined in Table Al of lLudwig et al.l dl999t) . By chang- 



ing the value of fa, the convective efficiency is modified in the 
transition between the optically thin and thick regions, which is 
the region of line formation. Therefore, this extra parameter can 
have a major effect on the predicted structures and spectra. 

Table 1 : Adopted parameterisation of MLT 



Param. 


Alt. name 


Value 


1/H P 


a 


0.8 


a 




1 


b 




2 


c 




16 


d 


r(i+/ 4 /r?) 


r(i + 2/T?) 



Notes. / is the mixing-length, H p the pressure scale height, T the con- 
vective efficiency and r e the optical thickness of a convective eddy. 



In Paper I, the LHD models were computed with optically 
thick convective cells (fa - 0) while in the present work, LHD 
models are allowed to have optically thin cells (fa = 2). As a 
result the LHD code is now based on an implementation of the 
MLT theory that is exactly the same as in commonly used white 
dwarf models. The effect on the 3D log g corrections is discussed 
in Sect. 4. We point out that 3D simulations are completely un- 
affected by the MLT parameterisation. 



2.3. Numerical precision 

Finally, we had to increase the numerical precision in order to 
provide the best match between the LHD and standard struc- 
tures. We found that at least 200 depth points must be u s ed with 
the LHD code in order to recover the Tremblav et al] (1201 lab 
structures. Furthermore, we find with both codes that the equa- 
tion of radiative transfer must be solved for at least five angles. 
We stress that our results in terms of the numerical precision of 
the LHD code can not easily be transferred to the CO s BOLD 
code, because of the completely different treatment of convec- 
tion. Instead, we show in Sect. 4.3 that the sensitivities of the 3D 
structures to different numerical parameters is very small and we 
conclude that our C0 5 BOLD code setup provides the required 
high precision. 



3. 3D model atmospheres 

We have computed a set of twelve 3D DA model atmospheres 
with the C0 5 BOLD codefl covering the range 13, 000 > T eS (K) 
> 6000 at loggr = 8.0. Table 2 gives a summary of the proper- 
ties of these simulations. We have used the new EOS and opac- 
ity tables that were discussed in Sect. 2. We relied upon four 
different opacity tables, with the opacity bins sorted based on 
a 6000, 8000, 10,000 and 12,000 K ID reference model atmo- 
sphere, respectively. The implementation of the boundary condi- 



1 We rely on the January 2012 version of the code which is similar 
to the version described in Frey tag et al.l d2012l) . We utilise some new 
features that we describe in this section. 



3 



Tremblay et al.: Pure-hydrogen 3D model atmospheres of cool white dwarfs 



tions is described in detail in lFrev tag et al. f 20 1 2l see Sect. 3.2) 
and CO 5 BOLD solar models rely on the same conditions. In 
brief, the lateral boundaries are periodic, and the top boundary 
is open to material flows and radiation. The bottom layer is open 
to convective flows (and radiation) in all but the three hottest 
simulations, and a zero total mass flux is enforced. We specify 
the entropy of the ascending material to obtain approximately 
the desired T e ff value (derived from the emergent stellar flux). 
The entropy of the deep layers is increasing monotonically with 
r e ff , hence there is a unique relation between the entropy at the 
bottom boundary and the effective temperature of a simulation. 
This control of the T e ff of a model (i.e. the temporal and spatial 
average of the emergent radiative flux) is indirect and r e jf is not 
an input parameter In contrast, the bottom layer is closed with 
imposed zero vertical velocities (but open to radiation) for the 
three hottest simulations, in which case we impose the radiative 
flux according to the diffusion approximation. 

We typically started our simulations by scaling approx- 
imately another model from the grid. As long as the pro- 
gram runs smoothly , the initial structure does not matter much 
dFrevtag et al.l2 012) for space- and time-averaged quantities. We 
first computed grey models for 10 seconds, and then switched 
on the non-grey radiative transfer for 10 seconds or more. More 
details about the time evolution of our simulations is given is 
Sect. 3.3. 

We adopt a grid of 150 x 150 x 150 points in the x, y and z 
directions, where z is used for the vertical direction and points 
towards the exterior of the star. Compared to the simulations 
computed in Paper I, the vertical number of points was in- 
creased from 100 to 150. The geometrical dimensions are given 
in Table 2. We fixed the bottom layer at log tr = 3 for all 
simulations, well below the photosphere^. One exception is the 
third hottest model (T e g ~ 12,000 K) for which we recom- 
puted a more extended simulation with a bottom boundary at 
log tr = 3.1. The bottom of the convective zone was too close 
to the simulation boundary at log tr = 3.0, which was arti- 
ficially damping the convective velocities by a small amount. 
Nevertheless, we find that our two simulations with an unequal 
vertical extent show very little differences except for the sub- 
photospheric convective velocities. In all models, the top bound- 
ary reaches a space- and time-averaged value of no more than 
log tr ~ -5. In Table 2 we also show the number of pressure 
scale heights between the photosphere and the bottom boundary 
that are covered by the simulations. It demonstrates that convec- 
tive eddies reaching the photosphere are unlikely to be impacted 
by boundary conditions. The grid spacing in the z direction is 
non-equidistant. The horizontal geometrical dimensions are fur- 
ther in discussed in Sect. 3.4. 

In terms of solving the hydrodynamical equations, the 
van Leer slope limiter was used in Paper I, but we have 
now switched to a less dissipative 2nd-order reconstruction 
method. Furthermore, regarding the time integration sch eme, we 
now a dopt the comer-transport upwind (CTU) method (Colella 
fl990l) . 

It is well known that the flows in stellar atmospheres are 
characterized by very large Reynolds numbers (Re > 10 10 ), and 
white dwarfs are not an exception. This implies that the flows are 
highly turbulent, and that the turbulent kinetic energy is then dis- 
sipated into heat at the Komolgorov microscale (d ~ // p Re -3 ^ 4 ), 



2 For clarity, we use round numbers in the discussion of the models 
in Table 2, but all calculations are done with the exact T e g values. 

3 We define the photosphere as the line-forming regions as opposed 
to the atmosphere which is the full simulation. 



i.e. on scales much smaller than the resolution of our grid. As a 
consequence, only the largest flow structures are resolved, and 
the small-scale kinetic energy is dissipated at the grid scale, ei- 
ther by the Roe solver itself or by an additional artificial tensor 
viscosity. The granules are the large flow structures that con- 
tribute to the global convective energy exchanges in the atmo- 
spheres, and as long as they are well resolved, the numerical 
treatment of the energy dissipation at smaller scales has very 
little effect on the predicted structures. Hence, viscosity is not 
a free parameter like, for instance, the mixing-length parame- 
ters. In Paper I, artificial viscosity was added, but with our new 
setup of the numerical schemes, we could compute all simula- 
tions without this explicit viscosity. This allows resolving a bit 
further the flows for which the spatial scale is close to the reso- 
lution of the box. Finally, we note that the radiation pressure is 
negligible for these models. 

Figs.[3]to|5]present the final snapshots for three of our sim- 
ulations at T eff ~ 8000, 10,000 and 12,000 K. The first series of 
plots (left panels) present the atmospheric structures in terms of 
temperature contours and velocity fields as a function of the geo- 
metrical depth (z) and one of the horizontal direction (x). We also 
show surfaces of constant Rosseland optical depth over which 
the mean structures are computed throughout this study. 

In all snapshots, significant velocities are observed in the up- 
per layers (tr < 10~ 3 ) even though these regions are stable to 
convection. The presence of oscillating waves is the main expla- 
nation for this behaviour, and it is difficult to differentiate these 
waves from the exponentially decaying convective flows (over- 
shoot layer), which transports energy outside of the convective 
zone (see Sect. 3.2). Another obvious observation is that the ve- 
locity field is significantly different for the hotter 12,000 K sim- 
ulation. It is seen that the downdrafts are concentrated in narrow 
lanes in comparison to the cooler models. Furthermore, roughly 
the lower third of the 12,000 K simulation, in terms of geomet- 
rical depth, is stable to convection. However, large downdrafts, 
such as the one on the far right of Fig. [5] reach high tr values. 
The root-mean-square (RMS) vertical velocities actually remain 
well above zero up to tr ~ 10 3 because of these downdrafts. 

The second series of plots (right panels) present the emer- 
gent frequency integrated (bolometric) intensity. It is shown that 
the granulation size is changing significantly over the range of 
Jeff- Once again, the 12,000 K simulation stands out with nar- 
rower cool downdrafts and smoother hot cells. The size, shape 
and contrast of the granulation patterns will be discussed in more 
details in Sect. 3.4. 

3. 1 . Mean structures 

We have shown in Sect. 2 that by using the proper EOS and 
opacity tables, it is possible to reproduce standard white dwarf 
model atmospheres with the LHD code. Given that CO s BOLD 
employs the same microphysics and radiative transfer schemes 
as LHD, we can conclude that CO s BOLD reaches a similar pre- 
cision, and possibly a better accuracy due the improved convec- 
tion treatment, in comparison to the current ID models. We now 
compare direc tly the results from th e 3D simulations with the 
ID models of lTremblav et al.l d201 lah . 

It was demonstrated in Paper I (Fig. 4) that in the case of 
DA white dwarfs, 3D spectral synthesis produces the same re- 
sult as computing one (3D) spectrum from a single (3D) struc- 
ture, which is the temporal and spatial average of a 3D simula- 
tion over surfaces of constant Rosseland optical depth. The mean 
temperature and pressure structures, and the atmospheric param- 
eters (r e ff and \ogg), are the quantities necessary to produce 
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log 9 


xX y X z 


ln(P bot/ Pphot)" 


Time 6 


3D log g corr. 


(K) 




(km) x (km) x (km) 




(s) 




5997 


8.0 


1.19x1.19x0.34 


5.0 


100 


0.00 


7012 


8.0 


1.19x1.19x0.38 


4.2 


60 


-0.07 


8032 


8.0 


1.40x1.40x0.43 


3.7 


60 


-0.16 


9035 


8.0 


1.58x1.58x0.51 


3.3 


10 


-0.26 


9520 


8.0 


2.06x2.06x0.62 


3.8 


10 


-0.25 


10018 


8.0 


2.10x2.10x0.76 


4.3 


10 


-0.24 


10530 


8.0 


2.24x2.24x0.88 


3.6 


10 


-0.22 


11004 


8.0 


2.68x2.68x1.19 


4.0 


10 


-0.15 


11531 


8.0 


3.65x3.65x1.69 


3.6 


10 


-0.11 


12022 


8.0 


7.45x7.45x4.06 


2.0 


10 


-0.07 


12505 


8.0 


7.45x7.45x4.08 


1.9 


10 


-0.09 


12999 


8.0 


7.45x7.45x4.41 


1.8 


10 


-0.09 



Phot is the pressure at the bottom layer and _P p hot the pressure at tr = 2/3. ' ' Total stellar time. For the three cooler models, the total time 
includes the initial 2D run as described in Sect. 3.3. 



1.2 




x (km) 




Fig. 3: Snapshot of the 3D white dwarf simulation at T e g ~ 8000 K and loggr = 8. Left: Temperature structure for a slice in the 
horizontal-vertical xz plane through a box with coordinates x,y,z (in km). The temperature is colour coded from 14,000 (red) to 
3000 K (blue). The arrows represent relative convective velocities, while thick lines correspond to contours of constant Rosseland 
optical depth, with values given in the figure. Right: Emergent bolometric intensity at the top of the horizontal xy plane. The RMS 
intensity contrast with respect to the mean intensity is 7.6%. The length of the bar in the top right is 10 times the pressure scale 
height at r R = 2/3. 



model spectra from the 3D simulations, which is the objective 
of Sect. 4. We derived mean T and P values from the average 
of T 4 and P over surfaces of constant tr. We selected 12 ran- 
dom snapshots in the last 2.5 seconds of the simulations to make 
the temporal average. We determined the T e ff values identified 
in Table 2 from the mean emergent flux of the 12 snapshots se- 
lected for the tr average. The RMS variation of r e ff with time is 
always much less than 1%. In contrast, \ogg is a fixed quantity 
in the hydrodynamical equations. 

To understand the differences between (3D) and ID models, 
it is useful to look at other mean quantities. Therefore, we com- 
puted the mean entropy and RMS velocities at constant optical 
depth, and different energy fluxes described below at constant 
geometrical depth. Fig. [6] illustrates that the characteristic depth 
of the formation of H6, taken as a typical line, is within the range 
-2 < log tr < 0. Hence this region of the photosphere is the 
most critical for the predicted spectra. 



In Figs. [7] and [8] we compare the (3D) and ID structures 
in terms of the temperature and e ntropy, respectively. Th e ID 
models are drawn from the grid of iTremblav et al.l (1201 lab with 
the ML2/a = 0.8 parameterisation of the MLT Compared to the 
temperature, the entropy values are more sensitive to differences 
between (3D) and ID structures. Furthermore, the presence of a 
negative or positive entropy gradient (as a function of tr) shows 
whether a layer is stable or unstable to convection, respectively. 
It is observed that in stability terms, the top of the convective 
zone, and even the bottom of the zone for the three hottest struc- 
tures, are at similar optical depths in (3D) and ID models. In 
the upper part of the atmospheres, the entropy gradient (in ab- 
solute value) and the temperatures are always much smaller in 
(3D) structures, which is the result of convective overshoot. This 
overshoot layer is able to cool the upper layers because of a weak 
radiative coupling (see Sect. 3.3). In comparison, the ID models 
are stable to convection and in radiative equilibrium. 
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Fig. 4: Similar to Fig.[3]but for the simulation at T e ff ~ 10, 000 K. Left: The temperature is colour coded from 17,000 (red) to 5000 K 
(blue). Right: The RMS intensity contrast is 14.4%. 




Fig. 5: Similar to Fig.[3]but for the simulation at T s $ ~ 12, 000 K. Left: The temperature is colour coded from 60,000 (red) to 7000 K 
(blue). Right: The RMS intensity contrast is 18.8%. 



In Fig.|9l we present the ratio of the convective flux to the to- 
tal energy flux as a function of optical depth. The convective flux 
is the sum of the enthalpy, gravitational energy and kinetic en- 
ergy fluxes. The shape of the convective zones as described by 
the convective flux is also similar between (3D) and ID mod- 
els. One reason is that the overshoot in the upper layers acts 
in optically thin regions, and transports very little flux. For the 
three hottest models, however, the convective zone is systemat- 
ically deeper for the (3D) models in part because an overshoot 
layer is present at the bottom of the convection zone where it 
is able to transport energy. This is in qualitative agreement with 
the results from pulsation studies in which a slightly more effi- 
cient MLT parameterisation (ML2/ g _ 1.0) is used to describ e 
the bottom of the convective zones (Fontaine & Brassard 2008). 
We note that larger than average differences in the (3D) and ID 
temperature structures at log tr > 1 for the hot models are due 
to the transition from convective to radiative layers. In general, 
the transition region from convective to radiative flux transport is 



smoother in (3D) structures according to Fig. which explains 
why (3D) structures have smoother temperature gradients in the 
photosphere as reported by Fig. |7] 

Interestingly, the (3D) and ID models are the most similar 
at the boundaries of our computed sequence and the maximum 
differences between ID and (3D) models in the photosphere are 
detected in the middle of our sequence at ~ 9000- 10, 000 K. 
At r e jf ~ 6000 K, convection is almost fully adiabatic every- 
where in the atmosphere, and since ID and 3D models are based 
on the same input microphysics, it is expected that the adiabatic 
gradient, hence the structures, will be the same. For the hottest 
model (T e ff ~ 13,000 K), the maximum convective flux becomes 
small, and the atmosphere is leading towards the fully radiative 
regime. Nevertheless, it is remarkable that the entropy gradients 
are still substantially different in the upper layers, which sug- 
gests that the convective overshoot is still significant at this tem- 
perature. However, part of the differences between the ID and 
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Fig. 6: Similar to Fig.[2]but in the wavelength range, with respect 
to the line center, of the red wing of H<5. 
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Fig. 7: Temperature structures versus log tr for (3D) (red, solid) 
and ID (black, dashed) model atmospheres. The temperature 
scale is correct for the 6000 K model (bottom curve), but other 
structures are shifted by 1 kK relative to each others for clarity. 



3D T e ff structures may also be caused by the distinct treatment 
of the radiative transfer. 

3.2. Convective velocities 

Another relevant aspect of the 3D simulations are the convective 
velocities. We present the RMS vertical and horizontal veloci- 
ties in Figs.[10]and[TT| Due to the symmetry of the simulations, x 
and y velocities are nearly identical and we only show the former 



Fig. 8: Entropy as a function of log tr for our sequence of (3D) 
(red, solid) and ID models (black, dashed), with the coolest 
model at the bottom, and the hottest model at the top. The en- 
tropy scale is correct for the 6000 K model, but other structures 
were shifted by 0. 1 unit for clarity. 



component. Our sequence of models exhibits smooth changes of 
the velocities between adjacent T e ff values, which further con- 
firms that our structures are well converged. The vertical veloci- 
ties in Fig.[l0]are the most interesting since we can relate them 
to the convective energy flux. The maximum vertical velocities 
are observed in the 1 1,500 and 12,000 K simulations and the val- 
ues decrease in a regular way in cooler and hotter models. The 
velocities in the coolest computation at 6000 K are only a small 
fraction of the maximum velocities observed in our sequence, 
which illustrates the fact that only mild convection is needed to 
transport the small total energy flux at 6000 K. 

The overshoot layer at the bottom of the convective zones is 
clearly seen in terms of the convective velocities. For the three 
hottest simulations, the velocities remain clearly above zero even 
at the bottom of the simulations (log tr ~ 3) while Fig.|8]shows, 
from the change of sign in the entropy gradient, that the layers 
become stable to convection at log tr ~ 1.1, 1.5 and 2.2 for the 
13,000, 12,500 and 12,000 K simulations, respectively. Clearly, 
the overshoot layers extend, in log tr, at least 1 dex and up to 2 
dex below the region unstable to convection. However, Fig. ^il- 
lustrates that the overshoot layers transport very little convective 
energy flux except in the vicinity of the convective zone. Clearly, 
the relevance of the overshooting depends on the physical pro- 
cesses that are studied. For instance, we believe that in diffusion 
or convective mixing studies, a prope r account of the over shoot 
layer is essential, as also concluded bv lFrevtag et al.l (ll996) from 
2D simulations. 

Figs. [TUl and [TTI would also suggest large overshoot layers 
above the convective zone, but we must be very cautious about 
this interpretation. Indeed, p-mode oscillations, with periods that 
are related to the geometry of the computational box (and bound- 
ary conditions), occur in all of our calculations. The effect of 
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Fig. 9: Ratio of the convective to the total energy flux as a func- 
tion of log tr for our sequence of (3D) (red, solid) and ID mod- 
els (black, dashed), with the coolest model at the bottom, and the 
hottest model at the top. The ratio is exact for the 6000 K model, 
but other structures were shifted by 0.5 up to 1 1,500 K, and by 
1 .0 above that temperature, for clarity. 




Fig. 10: Vertical RMS velocity v z as a function of log tr for our 
sequence of (3D) models from 6000 (lowest, red curve) over 
10,000 (top, light blue curve) to 13,000 K (dark blue curve). 



these oscillations is largely removed from all the mean quan- 



Fig. 11: Horizontal RMS velocity v x as a function of log tr for 
our sequence of (3D) models from 6000 (lowest, red curve) to 
13,000 K (dark blue). 



tities by the temporal average, except in the case of RMS ve- 
locities. The amplitude of the p-mode oscillations is larger in 
the upper layers, where the pressure restoring force is lower 
(b oc p~ 1 / 2 for undamped oscillations). In Fig. [10] the increase 
of the vertical velocities at log tr < -2 is caused by these oscil- 
lations. Methods hav e been proposed to filter these oscillations 
dLudwig et al.l l2002h in order to look at the exponential over- 
shoot more clearly. However, for the purpose of this work, the 
effect of overshooting into the upper layers can be more easily 
seen from the temperature and entropy structures. 

The turbulent pressure (pu?) is a significant fraction of the 
local gas pressure in convective white dwarfs. This fraction 
reaches values of ~10% for the hottest simulations of our grid. 
The turbulent pressure derived from the RHD simulations illus- 
trates the deviation from hydrostatic equilibrium. A proper ac- 
count of the turbulent pressure in hydrostatic ID models might 
help to improve the agreement between ID and (3D) structures. 

3.3. Characteristic time scales 

We note that three of our simulations (12,50 > r e ff (K) > 
1 1, 50 0) are within the ZZ Ceti instability strip dGianninas et al.l 
1201 ll) . It would be adequate to rely on our 3D structures as input 
for asteroseismic studies although the coolest ZZ Ceti star in our 
grid ( 1 1 ,500 K) is still convective at the bottom of the simulation. 
One option would be to compute RHD simulations that reach 
the bottom of the convective zone for all pulsating white dwarfs 
and then combine the RHD results with ID int erior structures. 
Such an application has been performed once by Gautschv et al.l 
(1996) and it would be appropriate to update this analysis with 
our improved and le ss dissipative 3D calculations. Alternatively, 
iLudwig et al.l ( 1 19991) have shown that it is generally possible to 
match a 3D structure including a convective bottom layer with 
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an unique ID structure model by relying on the asymptotic value 
of the spatially resolved entropy. 

The study of time-resolved white dwarf model atmospheres 
with hydrodynamical simulations is unparalleled by ID mod- 
els. We note that this time resolution is of interest for astero- 
seismic studies. It has been understood long ago that convec- 
tive processes were fairly rapid in ZZ Ceti white dwarfs com- 
pared to the typically observed pulsation periods in the range 
of 100-1000 seconds (Fontai ne & Brassard! I 2008). As a conse- 
quence, non-adiabatic pulsation codes have adopted the con- 
cept of instantaneous reaction o f the convection to the pulsa- 
tions (Bra ssard & Fontamal 19971) . The impact of time-resolved 
convection on asteroseismology of ZZ Ceti stars has recently 
been studied by Ivan Grootel et all d2012l) using perturbations in 
the mixing-length equations. However, they show that the differ- 
ences are small compared to the instantaneous reaction approxi- 
mation. 
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Fig. 12: Timescales for the (3D) structure of a 12,000 (top 
panel) and 8000 K (bottom panel) DA white dwarf as a func- 
tion of log tr. The different timescales are identified in the leg- 
end, and correspond to the advective timescale (black, solid), 
the Kevin-Helmholtz timescale (blue, dotted), the Brunt-Vaisala 
timescale (cyan, short dashed) and the radiative timescale (red, 
long dashed). 



In Fig. Q~2] we present four characteristic timescales as a 
function of the optical depth in the atmosphere of a 12,000 and 
8000 K white dwarf. The timescales have been evaluated un- 
der simplifying assumptions, and should therefore be taken as 
order-of-magnitude estimates only. The timescales are de fined 
as follow (see also lLudwig et al. 2002; Frevtag et al.ll2012l) : 



In 
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where St = — denotes the thermal expansion coefficient 

at constant pressure, V ac | the adiabatic gradient, V the actual tem- 
perature gradient, v c the RMS convective velocity, c p the specific 
heat per gram, <x the Stefan-Boltzmann constant, f is the char- 
acteristic optical depth of a disturbance of a size H p , with k, the 
mean bin opaci ty per gram, and w, the wei ght of the opacity bin 
(see eq. A.4 of lLudwig & Kucinskasll2012l) . 

First of all, the Brunt-Vaisala timescale (?bv) is related to the 
period of the oscillations due to the buoyancy force in regions 
stable to convection and the convective growth time for convec- 
tively unstable regions. The spikes observed for this quantity are 
related to the boundaries of the convective layers. The advective 
timescale (f a dv)> or the turnover timescale inside the convection 
zone, is related to the characteristic time for a convective mass 
element to move a distance of one pressure scale height. This 
timescale varies with depth from 10~ 2 to 1 second and it does 
not change significantly between the 8000 and 12,000 K simula- 
tions. 

The Kelvin-Helmholtz timescale (?kh) is the thermal energy 
content (per unit area) divided by the total energy flux and de- 
scribes the thermal relaxation (trough radiative diffusion). This 
timescale reaches values of the order of 15 seconds for the 
deeper layers of both the 8000 and 12,000 K models. However, 
in simulations below 12,000 K, the bottom layer is convective, 
and close to 100% of the flux is transported by convection. When 
the convective flux is much larger than the radiative flux at the 
bottom layer, the structure of the atmosphere is expected to ad- 
just on the advective timescale. For the three hottest simulations 
with a radiative bottom layer, we have verified that in practice, 
the total flux is stable after a few seconds of computation in all 
layers, perhaps due to the presence of convective overshoot and 
a correct initial guess of the structure of the deepest layers. 

The radiative timescale (f ra( j) refers to the characteristic time 
for the decay of local temperature perturbations through radia- 
tive transfer. The radiative timescale is of the order of the advec- 
tive timescale at 12,000 K, which is also the case for the 8000 K 
model, but only in the photosphere. The Peclet number in the 
photosphere, which is the ratio of the radiative timescale to the 
advective timescale, is changing from a value above to below 
unity between the cooler and hotter simulations, respectively. 
This implies that the evolution of convective cells in the photo- 
sphere will be largely influenced by radiation in the hot models 
in our grid, whereas for the cool models, convective cells are 
more directly governed by quasi-adiabatic expansion and com- 
pression. In the convective zone of the 12,000 K simulation, 
characteristic of a ZZ Ceti star, all timescales are lower than one 
second, which confirms that the reaction of the convective zone 
to changing conditions in the radiative layer just below will be 
rather instantaneous on the timescales of the pulsations. 

In the 8000 K model, the radiative timescale is rather large at 
the top and bottom layers. At the bottom layer, energy is trans- 
ported by convection and the large radiative timescale simply 
implies that the energy in the convective cells will hardly be lost 
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by radiation. The upper layers above log tr ~ -0.5 are, how- 
ever, stable to convection. Hence Fig.Q~2]suggests that these lay- 
ers may take up to 100 seconds to reach radiative equilibrium. 
The radiation field becomes weaker in cooler white dwarfs, and 
the decreasing line opacity creates upper layers that do not in- 
teract much with radiation. This contributes in creating the long 
observed radiative timescales . A similar phenomenon was o b- 
served for metal-poor dwarfs dAsplund & Garcia PerezlEOOlK 

To make sure that our simulations have relaxed in the upper 
layers, we performed non-grey 2D simulations for 50 seconds 
for the three models between 7000 and 9000 K, and for 90 sec- 
onds for the model at 6000 K. The upper layers never actually 
reach a radiative equilibrium like in the ID models. Instead the 
convective overshoot causes the entropy gradient in the upper 
layers to relax to a near-adiabatic structure, as seen in Fig. [8] We 
used the final 2D structures as initial conditions^ for the much 
more time consuming non-grey 3D simulations. 

All non-grey 3D simulations were run for 10 seconds, which 
ensures that we typically cover about 100 advective timescales 
in the photosphere. We have verified that all simulations are re- 
laxed in the last five seconds of computation and that they show 
no systematic and non-oscillatory change of their properties in 
all layers. We note that the four cooler simulations might still 
have an imprint of the 2D initial conditions in the upper layers 
because of the large radiative timescales. 

The numerical time step is limited by the global minimum 
of the radiative and Courant-Friedrichs-Lewy timescales (?cfl = 
Ax/[c s + v c ], where c s the adiabatic sound speed). The latter 
corresponds the travel time of the shortest wave across a grid 
cell of dimension Ax. Typically, our time steps are of the order 
10 - 10~ 5 second, and therefore the total number of steps is 
f sim /Af ~ 10 5 - 10 6 . 



teristic granulation is typically well resolved although because 
of the limitation of the Fourier analysis to identify granules and 
possible effects from the numerical parameters, the characteris- 
tic dimensions should be taken as estimates only. We note that 
we can not derive a single power law to describe the formation 
of smaller cells or sub-cells wi th a high wavenum ber, much like 
in main-sequence simulations (Ludwig et al. 200^). 
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Fig. 13: Mean power spectrum as a function of the horizontal 
wavenumber {In/A) averaged over 250 snapshots of the (bolo- 
metric) intensity map of the 10,000 K 3D simulation. 



3.4. Characteristic length scales 

We have previously discussed the vertical extent of our simula- 
tions which is naturally derived from our requirement to have 
the full atmospheres. However, the horizontal dimensions of the 
simulations are not well defined a priori. The 3D simulations 
have horizontal sizes, while such thing does not exist in ID 
plane-parallel models. From physical considerations in the for- 
mation, transport and dissipation of convective cells, one can as- 
sume that the vertical size of the cells will be of the order of 
the local pressure scale height. However, calculations of stellar 
RHD models have shown that granular patterns have larger hor- 
izontal dimensions o f the order of 10 times the local value of H p 
(iFrevtaget al.lll997l) . 

In Fig. [13] we present the power spectrum as a function of the 
horizontal wavenumber, for the emergent intensity in 250 differ- 
ent snapshots of the 10,000 K simulation (Fig.|4|gives one exam- 
ple of a snapshot). We display power per logarithmic wavenum- 
ber interva l for a more direct id entification of the power carry- 
ing scales (Ludwig et al. 2002). We remind the reader that the 
smallest wavenumber on the left in Fig. fT~3l represents a sinu- 
soidal pattern with one hot crest and one cool crest (compared 
to the mean intensity). Ideally, our simulations should include at 
least of the order of 3 x 3 hot cells to have a good representation 
of the atmosphere, and have a sufficient resolution for each cell. 
We only qualitatively ensured that it was the case since we can 
derive the power spectra only after the simulations. The charac- 

4 We copied all quantities of the 2D snapshots in the third dimension, 
except for the velocities which were taken from a previously computed 
3D snapshot. 



Fig. [T4lpresents the characteristic size of the convective cells 
in our sequence of simulations, derived as the wavelength of 
the peak of the power spectra which we found were well fit- 
ted by a parabola. It is clearly seen that the cell dimensions are 
not well scaled by the pressure scale height in the photosphere, 
also shown in the figure. The characteristic size of the granula- 
tion patterns increases much more rapidly with T e s than what we 
would expect from simple thermodynamics considerations. The 
convective velocities and the associated Mach number are signif- 
icantly increasing with T e e, and we suggest that it has an effect 
on the size of the convective cells by allowing for faster moving 
and thinner downdrafts. Furthermore, the timescales analysis of 
Sect. 3.3 suggests that the variation of the Peclet number is part 
of the explanation. The regime of regular bright and dark cells at 
cool temperatures (see, e.g., Fig. [3), where the advection domi- 
nates, change to the regime of large bright cells and narrow dark 
lanes at hotter temperatures (Fig. [5) when the radiation domi- 
nates. 

Also of interest is the intensity contrast of the granulation 
patterns with respect to the mean intensity. This quantity can be 
derived from the intensity snapshots such as the ones shown in 
Figs. [3][5] The RMS intensity contrast, a 3D effect that can not 
be easily predicted from existing ID structures, is presented in 
Fig. [15] It is seen that the contrast varies from 1 to 20% in our se- 
quence. In all but the three coolest models, the intensity contrast 
is in the 10-20% range. Very similar values of the intensity con- 
trast a re found for F, G and K dwarfs (see Fig. 15 of lFrevtag et al.l 
|20H . We find that the granulation is visually similar in the Sun 
and a ~1 1,000 K white dwarf. In cooler white dwarfs, the granu- 
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Fig. 14: Characteristic granular size (maximum of the mean 
power spectrum) as a function of r e ff for our sequence of 3D 
simulations. The points are connected for clarity. The solid line 
is eight times the pressure scale height at t r = 2/3. 



We mention that at 13,000 K, the intensity contrast is still 
high, while the convective flux is fairly low according to Fig. [9] 
It appears that because the total energy flux is large in the atmo- 
sphere, large temperature differences are necessary to transport 
even a small amount of convective flux. However, the differences 
between ID and (3D) models are small according to Fig. [7] This 
may not be entirely surprising, since the weak convective flux is 
unlikely to have a global impact on the largely radiative struc- 
tures. At the opposite cool end of our sequence, the 6000 K sim- 
ulation features an intensity contrast of about 1%. It suggests 
again that convection is very efficient to transport the small total 
flux in this regime with only very small temperature variations. 
In between these two extremes, the increasing total flux con- 
tributes to enhance the contrast. Furthermore, as is increas- 
ing, the Rosseland opacity is also increasing in the photospheres. 
This implies lower characteristic densities in the photospheres 
and higher convective velocities. This effect also contributes in 
enhancing the fluctuations. 

Finally, we recall here that despite the significant variations 
in the spatially resolved intensity or flux, results of Paper I have 
shown that the mean wavelength-dependent flux is almost ex- 
actly equal to the same quantity computed from mean (3D) 
structures. Hence the 3D fluctuations impact only indirectly the 
predicted spectra of white dwarfs through global modifications 
of the structures. In the following, the 3D simulations are there- 
fore represented by their mean structures. 



lation appears less structured in comparison to K main-sequence 
stars (see Fig. |3j. We note that cool white dwarfs have a wide 
entropy minimum, which implies that the upper boundary of the 
convective zone is not well defined and that granulation patterns 
may form in a more extended optical depth range than dwarfs. 
This is a possible explanation for the fuzziness of granulation in 
cool white dwarfs. A careful study of the granulation patterns as 
a function of T e ff, log g and metallicity for stars and white dwarfs 
would be useful to understand these structures. 
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Fig. 15: RMS intensity contrast divided by the mean intensity as 
a function of T e ff for our sequence of 3D simulations. The points 
are connected for clarity. 



4. Astrophysical applications and discussion 

4.1. Model spectra 

Our sequence of 3D model atmospheres can be used as input 
for spectral synthesis as done in Paper I. We have shown in 
Sect. 2 that our 3D simulations are based on the same input 
microphysics as the standard I D models, hence we ca n use the 
ID model atmosphere code of Tremblav et Id] (I201 lal) to com- 
pute (3D) spectra, from (3D) structures, that can be compared 
directly to observations. However, since our sequence of 3D cal- 
culations is limited to one log g value, we can not fit actual obser- 
vations and we will still rely on a (3D)-1D differential approach 
to look at the 3D effects. We compare in Fig.[l6]the (3D) and ID 
spectra for three characteristic simulations from our sequence of 
model atmospheres identified in Table 2. For clarity, only the 
blue wings of the lines are shown since the comparison is simi- 
lar for the red wings. We must conclude, like in Paper I, that the 
differences are fairly subtle between the predicted (3D) and ID 
spectra, and that we must be careful about drawing conclusions 
about the atmospheric parameter corrections. 

We find that the largest differences between (3D) and ID 
spectra are in the middle of our sequence, at T e ff ~ 10,000 K, 
which is also the region where the structures are the most dif- 
ferent (see Fig. |7]i and where the high-log g problem is the most 
significant (see Fig. 1 of Paper I). One surprising result is that 
the higher lines of the series are significantly impacted by 3D 
effects. Since the higher series members are formed in a nar- 
row region of the photosphere, one would expect that 3D ef- 
fects decrease for these lines. The reason for this behaviour is 
that the strength of the higher lines is sensitive to the non-idea l 
effects (iHummer & Mihalas|[l988l: iTremblav & Bergeronl l2009) 
in comparison to the lower lines. The non-ideal effects are in 
turn sensitive to the density in the atmospheres, and (3D) struc- 
tures have systematically cooler temperatures and higher den- 
sities in the upper layers due to the convective overshoot. These 
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Fig. 16: Comparison of the blue wing of six Balmer line profiles (Ho- to H8) calculated from (3D) structures (red, solid) and standard 
ID structures (black, dashed) for three models from our sequence identified in Table 2. T e ff values are shown on the different panels. 
All line profiles were normalized to a unity continuum at a fixed distance from the line center. All spectra were convolved with a 
Gaussian profile with a resolution of 6 A to represent typical observations. 



higher densities enhance the non-ideal effects and reduce the line 
strengths in comparison to the ID predictions. 

The cores of Ha and to a lesser degree H/3 are significantly 
deeper in (3D) spectra. These cores are formed very high in 
the atmospheres, where the (3D) structures deviate significantly 
from their ID counterparts due to the overshoot cooling. Since 
spectroscopic analyses of white dwarfs with ID models provide 
good fits to the line centers, the predictions from the (3D) spectra 
should be regarded with caution. We believe that the lower tem- 
peratures predicted by our 3D simulations are real, although we 
can not rule out that shortcomings in radiative transfer or miss- 
ing physics (e.g. magnetic fields or shock formation) might have 
an impact on this issue. We have verified with the TLUSTY code 
dHubenv&Lanzll 19951) that the ID NLTE effects are restricted 
to the very center of the lower lines (A/1 < 1 A) in this regime of 
r e ff and that they contribute to further enhance the absorption. 
Therefore, it seems unlikely that NLTE effects can help in ex- 
plaining the differences between (3D) and ID line cores. As a 
conservative measure, we remove the line centers (|AA| < 1.5 A) 
in our derivation of the 3D atmospheric parameter corrections. 

4.2. Application to the high-log g problem 

We derive here 3D atmospheric parameter corrections that are 
defined as the differences in T e g and log g when we fit the nor- 
malized line profiles of the (3D) spectra with our grid of stan- 
dard ID models. More specifically, the corrections are defined as 
1D-(3D), since we are interested in how much the atmospheric 
parameters of real stars would change when using (3D) model 
spectra. In Fig.fTTIand in Table 2, we present the 3D log g correc- 
tions. As in Paper I, we find that the corrections are negative for 
all spectra (i.e., the 3D simulations predict lower surface grav- 
ities), which is in the appropriate direction to correct the high- 



log g problem. We note that the 3D log g corrections are smaller 
in absolute value than those found in Paper I in the range of T e s 
in common between both studies, i.e. the four hottest simula- 
tions. This is a direct consequence of the fact that the ID LHD 
models used in Paper I had optically thick convective cells ev- 
erywhere in the atmosphere and were inconsistent with standard 
ID models (see Sect. 2.2). 

Also given in Fig. [17] are the observed shifts in log g de- 
rived from the spectroscopic a nalyses of the DA whit e dwarfs in 
the Sloan Digit al Sky Survey (Trembla y et al.ll201 lal) and White 
Dwarf Catalog (iGianninas et al.ll201 lb . These shifts correspond 
to the corrections required, in a bin of 1000 K around the simu- 
lation temperature, to match the mean mass value obtained from 
hot DA stars. We can see that the 3D corrections describe well 
the shape of the high-log g problem, with maximum corrections 
at r e ff ~ 10, 000 K where the problem is the largest. Also, the 3D 
corrections have about the right amplitude everywhere to solve 
the high-log g problem. Knowledge of the 3D log g corrections 
at loggr = 7.5 and 8.5, and ultimately the spectroscopic fit of 
observations will be necessary to further constrain the log g val- 
ues predicted by the 3D simulations. However, our results con- 
firm the conclusion of Paper I that ID MLT convection is the 
main reason for the high-log g problem, and that 3D model at- 
mospheres provide a more stable surface gravity distribution as 
a function of the effective temperature. 



We find that the 3D T e fi corrections are mild, with an aver- 
age of -230 K for the simulations in our sequence, hence the 3D 
effects are mostly log g effects. We believe that an actual com- 
parison with observations will be necessary to interpret further 
the r e ff shifts. 
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4.3. Sensitivity to numerics 

In order to understand the precision of our 3D simulations and 
the uncertainty of our log g corrections, we computed a series 
of 9 simulations with the same input parameters, except for one 
modification as given Table 3. All simulations were computed 
at r e ff ~ 10,000 K and \ogg - 8 for 5 seconds (at least one 
order of magnitude larger than both f ra( j and f a d v in the photo- 
sphere), with the converged model from our regular sequence as 
the starting model. We derived the mean temperature and pres- 
sure structures from 6 snapshots, and computed spectra using the 
same approach as we did for our regular sequence. In Table 3, we 
present how the 3D log g corrections are changed by our modifi- 
cations. 



o.io 

0.00 



-0.20 
-0.30 
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III 
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3D log g corrections 
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Observed in SDSS 
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Observed in WD Catalog 





















T e „ (K) 

Fig. 17: 3D \ogg corrections (red filled circles) as a function of 
r e ff. In comparison, we show the shifts in surface gravity re- 
quired to obtain a stable log g distribution as a function of r e ff 
using the samples of the SDSS (open circles) and Gianninas et 
al. (201 1; open squares). The points are connected for clarity. 



rived in this paper, these could be valid alternative input param- 
eters for models in the regular sequence. We have A) increased 
the horizontal geometrical dimensions x and y by 25% (hence 
lowered the horizontal resolution), B) decreased the vertical res- 
olution by using 100 grid points instead of 150, C) cut a section 
of atmosphere by changing the maximum value of log tr from 
3.0 to 2.8, D) changed the top boundary condition by increas- 
ing by 10% the temperature of the incoming flow, E) used the 
Van Leer slope limiter, as in Paper I, instead of the less dissipa- 
tive 2nd order method, F ) used the mo re aggressive Piecewise 
Parabolic reconstruction (IColellalll984l) and G) added artificial 
viscosity as in Paper I. According to Table 3, all of these mod- 
ifications have no significant effect on the predicted gravities. 
This is a very important result, which confirms that the unre- 
solved turbulent energy dissipation is not an issue for the pre- 
dicted Balmer lines. 

The simulations H and I include modifications that presum- 
ably result in a significantly less accurate account of the physical 
processes in comparison to the models in our regular sequence. 
We computed a 2D simulation by removing one of the horizontal 
direction. We find that the log g correction is changed by +0.06 
dex, which is a moderate but significant effect. This result is 
meaningful since 2D models have been used in this work to pro- 
vide relaxed upper layers for cool simulations (see Sect. 3.3). 
In fact, since 2D simulations are an order of magnitude faster to 
compute than 3D models, they are an interesting resource for ap- 
plications where the geometrical dimensions or calculation times 
are larger than those in our sequence of 3D simulations. While 
these 2D calculations may not be as accurate as the 3D setup 
to predict the Balmer line spectra, they still appear to provide a 
very good physical description of the atmospheres. We note that 
early works on RHD models of white dwarfs were done with 2D 
simulations. Our result shows that their 2D approximation was 
adequate to derive white dwarf properties. 

We close this section with the sim ulation I where we r elied 
on the opacity table first computed bv lLudwig et al.l d 19941) and 
also used in Paper I. We find that the log g correction is only 
changed by +0.03 dex. Therefore, we conclude that the opac- 
ity binning procedure was already very well tuned in previous 
works on RHD models of white dwarfs. 



Table 3: Sequence of 3D simulations with alternative parameters 



Description 


log g shift 


A) Surface area A XIJ = 1.56 A^o 


0.00 


B) Vertical grid points n, = 2/3 n ; ,o 


0.00 


C) logT,„ flv = 0.93 logTWa.o 


-0.01 


D) Top layer inflow T = 1.1 x To 


0.00 


E) Van Leer slope limiter 


-0.00 


F) Piecewise Parabolic reconstruction 


-0.02 


G) Artificial viscosity 


0.00 


H) 2D calculation 


0.06 


I) Original opacity table (Ludwig et al. 1994) 


0.03 



Notes. All models with T eS ~ 10, 000 K and loggr = 8. The subscript 
refers to the standard model parameters. 

The first category of modifications (simulations A to G) can 
be classified as tolerable, i.e. from the physical constraints de- 



5. Conclusion 

We computed a sequence of twelve 3D model atmospheres with 
the CO 5 BOLD radiation-hydrodynamics code. These simula- 
tions were made in the range 13,000 > T e s (K) > 6000 and at 
log g = 8 for pure-hydrogen atmospheres. We relied on EOS and 
opacity tables that have the same microphysics as in standard ID 
models of white dwarfs, and we have verified that our models are 
expected to have the same precision as prevalent ID structures. 
As a consequence, our sequence of 3D simulations is now ex- 
pected to predict realistic absolute properties and can ultimately 
be compared with observations. We have demonstrated that the 
3D simulations depend only weakly on numerical parameters by 
running a set of alternative simulations with modified parame- 
ters. This is unlike ID models for which the free parameters in 
the mixing-length theory have a significant impact on the pre- 
dicted structures. 

We derived (3D) spectra of the Balmer lines that we com- 
pared with those predicted from standard ID white dwarf mod- 
els. We found that the 3D \ogg corrections are the largest at 
around T e « ~ 10,000 K and that they have the right ampli- 
tude as a function of r e ff to draw the conclusion that the source 
of the long-standing high-log g problem is the inability of the 
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mixing-length theory to properly account for the convective en- 
ergy transport. We will follow this study with the calculation of 
sequences of simulations at \ogg = 7.0, 7.5, 8.5 and 9.0. It will 
then be possible to apply the 3D model atmospheres in the spec- 
troscopic analysis of DA white dwarf samples. 
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